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Abstract 


High-frequency sampled multivariate continuous time autoregressive moving average pro¬ 
cesses are investigated. We obtain asymptotic expansion for the spectral density of the sampled 
MCARMA process (V„A)nez as A f 0, where (It)tgR is an MCARMA process. We show that 
the properly filtered process is a vector moving average process, and determine the asymptotic 
moving average representation of it, thus generalizing the results by Brockwell et al. [9, 10] in 
the univariate case to the multivariate model. The determination of the moving average repre¬ 
sentation of the filtered process, important for the analysis of high-frequency data, is difficult 
for any fixed positive A. However, the results established here provide a useful and insightful 
approximation when A is very small. 

Keywords: Multivariate continuous time autoregressive moving average (CARMA) process; 
spectral density; high-frequency sampling; discretely sampled process. 

AMS 2000 Subject Classification: Primary 60G12, 62M10, 62M15; 


Secondary 60G10, 60G25, 60G35. 


1 Introduction 

The main object of this paper is the multivariate continuous time autoregressive moving average 
process (MCARMA) in d dimension, which we define as follows. Let 


P{z) = IdzP + AizP-^ + ... + Ap, 
Qiz) = Boz'i + Biz'i-^ + ... + Bg 


be the autoregressive and moving average polynomial respectively, Ai,Bj € M^, and Id is the 
d-dimensional identity matrix. The set of m x n real matrices is denoted by and for 

m = n. The driving process is a two sided d-dimensional Levy process {Lt)t € M, that is 



Li(t), t > 0, 

-L2{-t-), t < 0, 


where Li{t), L 2 {t), t > 0, are independent identically distributed (one-sided) d-dimensional Levy 
processes, such that ELi(l) = 0, and E||Li(l)|p < oo, with || • || being the usual Euclidean norm. 

‘This research was funded by a postdoctoral fellowship of the Alexander von Humboldt Foundation. 


1 



The covariance matrix of Ti(l) is For definition and properties of Levy processes we refer to 
Bertoin [3]. 

The continuous-time analog of discrete time ARMA equations is the differential equation 

P{D)Yt = Q{D)DLt, 


with D being the differential operator with respect to t. Since Levy processes are not differentiable, 
this is meant as the following state space representation. 

The d-dimensional stochastic process Y is an MCARMA process with autoregressive and mov¬ 
ing average polynomial P and Q respectively, if 


dG{t) = AG{t)dt + BdLt, 
Yt = CG{t), teM, 


where 


A = 


/ 0 h 0 

0 0 Id 

0 0 0 

\, Ap Ap—i Ap—2 


0 \ 
0 

Id 
-All 


and 


G Mpdi C — {Id, 0,0,... , 0) G Md^pd, 


B={l3l,l3^,...,/3;)' GMpd,d, I3p-j = 


1 Aifip—j—i Bq—j, 0 Y j < q, 
0, j > q- 


( 1 ) 


Let Ai,..., Xpd denote the eigenvalues of A, which is the same as the of zeros of detP( 2 ;), see Lemma 
3.8 in Marquardt and Stelzer [19]. It is well-known (see Brockwell [5]) that in one dimension a 
stationary causal solution exists if and only if the zeros of detP( 2 ;) have negative real parts. Under 
this condition strictly stationary causal solution of the MCARMA equation (1) exists, see Schlemm 
and Stelzer [21, Proposition 3.2]. Therefore throughout the paper we assume that 


the zeros of detP{z) have negative real part. 


Under this assumption the process Y can be represented as a continuous time moving average 
process ([19] (3.38)~(3.39)) 

Ft = / g{t - s)dLs, t G M, (2) 

J —OO 

where 

g{t) = — P{\x)~^Q{\x)dx. 

J — CO 

By Lemma 3.24 in [19], the assumptions on the eigenvalues of A implies that g vanishes on the 
negative half-line, that is our process is causal. In this paper we use representation (2). Since we 
are only interested in second-order properties of the process Y, the integral in (2) is understood 
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in the L^-sense. For the same reason our results remain valid in a more general setup, when the 
process {Lt)t£R has stationary orthogonal increments. 

These processes were introduced by Marquardt and Stelzer [19] in 2007 as a continuous time 
analog of vector-valued ARM A processes, and multivariate extension of continuous time ARM A 
models. The state space representation of these models were investigated by Schlemm and Stelzer 
[21], in particular the definition given above is from [21]. Gaussian CARMA processes date back to 
Doob [14] in 1944, while more general Levy-driven CARMA models were introduced by Brockwell 
[6] in 2001. These models are important tools for stochastic modeling and have a wide range of 
applications, e.g. in financial mathematics, electricity markets, and in turbulence. Andresen et 
al. [1] used Gaussian CARMA processes to model the short and forward interest rate. It is worth 
to mention that the Vasicek model corresponds to the Ornstein-Uhlenbeck process, which is the 
CARMA(1,0) process. Todorov and Tauchen [22] applied Levy-driven CARMA processes to model 
stochastic volatility in finance. For a review on applications of Levy-driven time series models in 
finance we refer to Brockwell and Lindner [4]. For spot prices in electricity markets Benth et al. [2] 
proposed a model with seasonality, consisting of a deterministic seasonality, a CARMA factor 
driven by a non-Gaussian stable Levy process, and a non-stationary longterm factor given by a 
Levy process. Brockwell et al. [10] applied CARMA processes to model high-frequency sampled 
turbulence data. 

Investigations of MCARMA models have become more active in the recent years. Ergodic 
and mixing properties of equidistantly sampled MCARMA processes were studied by Schlemm 
and Stelzer [21]. Fasen [15] investigated asymptotic properties of high-frequency sampled models, 
and gave parameter estimation, while Fasen and Kimmig [16] studied information criteria for 
MCARMA models. Recovery of the driving Levy process was treated by Brockwell and Schlemm 
[12] and by Ferrazzano and Fuchs [17]. 

In the present paper we consider high-frequency sampling of an MCARMA process and inves¬ 
tigate the characteristics of the resulting process = YnA- As A 4- 0 we obtain the asymptotic 
moving average representation of the filtered process, thus extending the results in Brockwell et 
al. [9], and partly in [10] to the multivariate setup. Our results slightly generalize the moving 
average decomposition of Schlemm and Stelzer [21, Theorem 4.2], and gives further information 
on the decomposition as A 4. 0. The determination of the moving average representation of the 
filtered process, important for the analysis of high-frequency data, is extremely difficult for any 
fixed positive A even in the one-dimensional case. However, the results established here provide a 
useful and insightful approximation when A is very small. 

In the next section we give two representations of the spectral density of the sampled process. 
The first one in Proposition 1 is a Taylor expansion in A, which we use to prove the asymptotic 
moving average representation in Theorem 1. The second one in Proposition 2 and 3 allows us 
to show that the filtered process is a moving average process. Section 3 contains the main result, 
the moving average representation of the filtered process. In Section 4 the simplest non-univariate 
example is spelled out in detail, showing the difficulties to obtain higher-order approximations as 
is possible in the one-dimensional case. All the proofs are gathered in the last section. 
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2 Spectral density of the sampled process 


It is well-known that MCARMA processes have spectral density function, which is 


1 


/y(A) = —P(iA)-^g(iA)SiQ(iA)*(P(iA) 


— 1 \* 


A E 


where A* is the Hermite-transpose of the complex matrix A; see (3.43) in [19]. Put 

R{z) = P{z)-^Q{z)ElQ{-z)'^{P{- z)-^y, zgC. 


( 3 ) 


(4) 


Note that the components of the matrix P{z)~^Q{z) are meromorphic functions, which can have 
poles only at the zeros of detP( 2 ;). Due to our assumptions the poles have negative real part. 

We frequently use the following simple facts about residues. If h{z) is a meromorphic function 
and p is a closed curve, which encircles the poles of h, then its residue at infinity is defined as 

Kes{h{z),oo) = -r [ h{z)dz, 

Jp 

moreover, it can be computed as 

Res(/i(z),oo) = —Res(z“^/i(z“^), 0). 

For a matrix M (z) with rational function entries let deg M{z)ij denote the degree of the numerator 
minus the degree of the denominator in the (i,j)-th element of M. Finally, for a square matrix A 
its adjugate is the transpose of its cofactor matrix, that is adjA = C~^, where Cij = Aij, 

and Aij is the determinant of the matrix that results from deleting row i and column j of A. 
When A is non-singular then A~^ = adjA/detA. 

From now on, let us fix a closed curve p C (—oo,0) x iM in the left half complex plane, which 
contains the poles of R{z) (the zeros of detP( 2 ;)). For any non-negative integer k introduce the 
notation 

f z^R{z)dz = —2TriIles(^z^R{z),oo^=:2TTiQk- ( 5 ) 


IpU-p 


Since 


P '(^) = have 

Riz) = 


1 


-adiP{z)Q{z)J:L{adiP{-z)Q{-z)) 


T 


detP( 2 ;)detP(—z) 

Note that deg (adjP(z)) < {d — l)p and degg(z) < q, thus we may write 

{d-l)p+q 

(adjP(z))g(z) = ^ Sjzy 

j=0 


( 6 ) 


where Sj G M^, j = 0,1,..., {d—l)p+q. Since the polynomials on the main diagonal of adjP(z) are 
of degree {d— l)p, otherwise the degrees are strictly less than {d — l)p, we obtain that = 

Bq. With this notation 


2[(d-l)p+q] ^ 

5(z) = adjP(z)g(z)SL(adjP(-z)g(-z))^ = ^ SjP, 

j=0 


(7) 
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where, the first coefficients are 


From ( 6 ) we see that for j = 0,1,... 

Sjj = 5 * 2 j and = ”* 527+1 • 

In the latter matrices the main diagonal is 0. In particular, they are 0 in the one-dimensional case. 

From (7) we see that each component of the matrix z^R{z) decreases at least as z~‘^ for k < 2{p— 
q) — 2, thus from the definition of the residue at infinity follows 0^ = 0. For k = 2{p—q) — l it is easy 
to determine the limit matrix from the definition of R{z), which is 02 (p_q)_i = 

To determine further coefficients note that 


OO 




z 


( 8 ) 


As detP(z) = — Aj) we obtain 


z-^R{z-^) = 


(^ — \Yd^2pd-l 


2[{d-l)p+q] 

E 

7=0 


S^z-K 


This formula implies a linear recursion for the coefficients 0^, which in a special case is spelled 
out in Section 4. 

Let us define the coefficients Ckipj) for a; 7 ^ 0, a; € (— 7 r, 7 r), via the series expansion 

1 °° 

1 = E kl < l^^l- (9) 

k=0 

Let /a denote the spectral density matrix of the sampled process = Y{nA). In the following 
we obtain a Taylor-expansion in A for the spectral density matrix. 

Proposition 1. ITe have for uj ^ 0, u G (— 7 r, 7 r) 


/a(w) = E (-^)^0fcCfc(w). 

^ k=2{p-q)-l 


From the leading term, we obtain 

Corollary 1. As A J, 0 

/A(a;) = (-If-''- ^ >’ c2ip-,)-i{u;)BoJ:LBl 

In particular, the first order approximation is real. 
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By Lemma 1 below C 2 {p-q)-i{ijj) is real, which implies that the first order approximation is 
real. 

We also give another representation of the spectral density, from which the moving average 
representation (17) below follows. 


Proposition 2. Assume that Ai,... ,\pd are different zeros of detP{z), i.e. each has multiplicity 
one. Then 


fA{i^) 


1 e 
dvr 

e=i 


0 ({£)~^) — e^^‘a{£) + e 
cosh AA^ — cos uj 


( 10 ) 


where the coefficient matrices a{£) G Md come from the partial fraction decomposition of R{z), i.e. 


R{z) 


S{z) 

detP( 2 ;)detP(—z) 


V («(^) , m \ 

-z-Xe) 


( 11 ) 


The further assumption on the multiplicity of the zeros is not necessary but it makes the 
formulas simpler. The following proposition gives the spectral density in the general case. Note 
that, as an abuse of notation, now Ai,..., A^ are the different zeros of detP( 2 :). 

Proposition 3. Let Xi,..., Xm be the different zeros ofdetP{z) with multiplicity izi,... ,1/^. Then 


1 e ^‘^{a{£,l) — a{£,Vff) — e^^^a{£,l) + e 


£=1 

m A/ 

j 


cosh AXf — cos cj 


^^^s^,Xi)a{£,jy + .sf{-uj,Xe)a{£,j) 
e=i j=2 


( 12 ) 


(cosh AAf — cos cj)! 


where s^{(jj,Xe) are trigonometric polynomials of uj of degree j — 1, whose coefficients depend on 
A and Xi, and the coefficient matrices a{£,j) € Md eome from the partial fraction decomposition 
of R{z), i.e. 


s{z) ^ V V 

detP( 2 ;)detP(- 2 ;) - Xf)i {-z - Xf)i J ’ 


(13) 


3 Moving average representation 

Recall that Ai,..., Xpd denote the zeros of detP( 2 ;). Define the polynomial 

pd 

^^(z) = (l - , 

i=l 

and consider the filtered process 

= d>^(R)y„^, (14) 
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where B is the backshift operator. Note that, whenever A* is complex, its complex conjugate is 
also a root of detP( 2 ;), thus the polynomial has real coefficients. The power transfer function 
(see Theorem 4.4.1 in Brockwell and Davis [7]) of the filter is 



pd 

(^1 - 

i=l 


pd 

(cosh AXj — cos uj). 
i=i 


(15) 


When the zeros of detP(z) have multiplicity one, from (10) and (15) we see that the spectral 
density of the filtered process X‘^ is 


. , , -rr. c A ^ ^ ^ e-^^{a{e) - a(£)^) - e^^^a{£) + 


fmAi^) = 


dvr 


pd _i, 

]^(cosh AXj — cos uj) ■ 
j=i e=i 


cosh AA^ — coscu 


2Pdp^ r 

- TT (cosh AAj — cosoj) e~'‘^{a{i) — a{lY) — 


dvr 


£=1 


(16) 


This is clearly a trigonometric polynomial of degree less than or equal to pd. However, the coeffi¬ 
cient of (cosw)^*^ and of (cosctj)P'^ ^ sinw is a multiple of 

pd 

- a(ef). 

e=i 


Since in the partional fraction decomposition in (11) deg5(2;) < 2{{d — l)p -|- q), the sum above 
is necessarily 0. Therefore /maC*^) i® ^ trigonometric polynomial of degree strictly less than pd. 
In the general case, when the zeros are not necessarily different, by (12) we obtain a similar 
representation. Thus, we have shown the following generalization of Theorem 4.2 by Schlemm and 
Stelzer [21]. 

Corollary 2. For any A > 0 the filtered process X‘^ is a d-dimensional moving average process 
of order less than pd; i.e., there exist a matrix polynomial with degree less than pd and a white 
noise sequence Z^, such that 


x^ = <k^(H)y,f = (17) 

Schlemm and Stelzer also showed that the matrix polynomial can be chosen such that 
det'I''^(z) has no zeros in the closed unit disc. Note that in Theorem 4.2 in [21] it is assumed that 
the zeros of detP( 2 ;) are different, while the corollary above holds in general. In the one-dimensional 
case this result was shown by Brockwell and Lindner [11, Lemma 2.1], see also Proposition 3 by 
Brockwell et al. [8]. Our aim in this paper is to obtain further information on the sequence as 
A ^ 0. 

In order to state the asymptotic result for the moving average process we need a lemma about 
the coefficients Ck{u}). 
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Lemma 1. There exists polynomials qk-i,rk-i of degree k — 1 with real eoeffieients sueh that 


{2k - 1)! [2(1 - cosa;)]''c 2 fc-i(a;) = (-1)^ g'fc-i(cos w), 

(2/c)! [2(1 — coscj)]^’''^ C 2 fc(a;) = (—1)^ isinw • rfc_i(cos w). 

Moreover, 

k—1 k—1 

qk-i{x) = {-l)^-^2^~^ JJ(1 - X - C 2 k-i,j), and f 2 k-i,j = (2A: - 1)!2“(''"^), 

1=1 1=1 

and 

k—1 k—1 

rfc_i(x) = (-1)''"^2*'JJ(1 - X - ^ 2 fcj), and Y{^ 2 k,j = {‘^k)l2~^ . 

1=1 1=1 

For the zeros of these polynomials ^ 2 k-i,j,C 2 k,j 0 (0; 2), j = 1, 2,..., A: — 1. 

The first few polynomials and the numerical values of the corresponding roots are 

Qoix) = 1 , 

gi(x) = 2(x + 2),{-2}, 

q 2 {x) = 4(x 2 + 13x + 16), {-11.623, -1.377}, 

qsix) = 8{x^ + 60x2 297 x + 272), (-54.657, -4.141, -1.202}, 

q^{x) = 16(x^ + 251x^ + 3651x2 + 10841x + 7936), (-235.705, -11.59, -2.579, -1.126}, 
g 5 (x) = 32(x® + 1018x^ + 38158x^ + 274418x2 + 580013x + 353792), 

(-979.322, -30.003, -5.615, -1.973, -1.087}, 


ro(x) = 2, 

ri(x) = 4(x + 5), {-5}, 

r 2 (x) = 8 (x 2 + 28x + 61), (-25.619, -2.381}, 

r 3 (x) = 16(x^ + 123x2 + lOllx + 1385), (-114.258, -7.014, -1.728}, 

r 4 (x) = 32(x^ + 506x3 11706x2 + 50666x + 50521), (-481.928, -18.784, -3.832, -1.457}, 

r 5 (x) = 64(x® + 2041x^ + 118546x3 + 1212146x2 + 3448901x + 2702765), 

(-1981.48, -47.391, -8.116, -2.697, -1.315}. 

One sees that the polynomials have real roots, moreover, the roots have the interlacing property. 
However, we cannot prove this. Since the zeros of an orthogonal polynomial sequence have this 
property (see Chihara [13, Theorem 1.5.3]) it is tempting to think that the polynomial sequences 
above are orthogonal. However, it is easy to check that the recursion in Favard’s Theorem ([13, 
Theorem 1.4.4]) does not hold even for the first terms, therefore, neither of the two sequences of 
polynomials is orthogonal with any weight function. 

For G C let us define ??(^) = 1 — ^ ± ~ where the sign is chosen so that \ri{f)\ < 1. 

Now we can state the main result of the paper. 



Theorem 1. The moving average process has the asymptotie form 

p-q-l 

= (/, - n (1 - v{( 2 ip-g)-i,j)B)Zn, Zn ~ WN(0, Sz), as A f 0, 

i=i 

where 

A2{p-q)-l 

Sz =-i- . 

{2{p-q)-l)lUU-^m2ip-q)-ij)\ 

We note that r]{f 2 {p-q)-i,j) niight be non-real (although we conjecture that they are all real 
valued), in which case r]{f 2 {p-q)-i,j) also appears in the product, which means that the coefficients 
in the moving average expansion are real. 

It is interesting to observe that up to the first order asymptotic the matrix moving average 
polynomial is in fact a scalar polynomial, and the covariance structure only appears in the covari¬ 
ance matrix Sz of the white noise. Thus Theorem 1 has the same form as the first order version of 
Theorem 1 in [10]; see also Theorem 1 in [9]. Finally, we mention that the corresponding higher- 
order version of Theorem 1, the analog of Theorem 1 in [9], seems hopeless to prove. The proof 
breaks down on the factorization of the spectral density of the moving average process, since for 
matrix spectral density no factorization holds in general; compare Theorem 10 and 10’ in Hannan 
[18, Chapter II]. 


4 An example 

Let us consider the simplest possible non-univariate case. That is p = 1, q = 0, d = 2. Then 

P{z) = IdZ + Ai, Q{z) = Bq, Ai,Bo£ M 2 . 

Moreover, 

adiP{z)Q{z) = zBq + adJHiHo, 

that is in formula (6) Si = Bq,So = adjHiHo. Furthermore, in (7) we have 

1 


R{z) = 


detP(z)detP(—z] 

1 

detP( 2 ;)detP(—z] 


■izSi + So)J^L{-zSl +S^) 
■ (^ 2^2 + + 5o) , 


(19) 


with 


52 = -BqJ^lB^ , 5i = (adJHi)' - adj^iHoS^Ho' , 5o = adjHiHoSiHo' (adJHi)' . (20) 

Assume that the zeros of detP(z) are different. From (19) we can compute the matrices in the 
partional fraction decomposition. Standard calculation gives that the matrices in Proposition 2 
are 

1 


a(l) = 
a{2) = 


2 Ai(Af — A 2 ) 

-1 

2A2(A? - Ai) 


(52Af + 5iAi + 5o), 
(52Ai + 5iA2 + 5o). 
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Then using formula ( 16 ), lengthy but straightforward calculation gives 

sinhAiA sinhA 2 A\ ~ . , , » ^ , a 

---j + 52 (Ai sinh AiA - A 2 sinh A 2 A) j 

+ isino; • 5 i(cosh AiA — cosh A2A) 

cosh Ai A • sinh A2A cosh A2A • sinh AiA 

A 2 Ai 

+ S2{X2 cosh AiA • sinh A2A — Ai cosh A2A • sinh AiA) 



/ma(w) = 


2e^(^i+^2) 

27r(Af - A^) 


COSCJ 




The corresponding process is MA(1), and according to Theorem 10’ of Hannan [18] there is a 
positive symmetric real matrix Tq and a real matrix Ti, such that 

/ma(^) = + 'hie‘‘")('ho + 


After short calculation one sees that the first order expansion is ~ {I — B)Zn, with covariance 
matrix = ABqTjlBq , as we have shown in Theorem 1. However, in general determining exactly 
the matrices is difficult. 

We can also use Proposition 1. Combining (19) and (8) the @k matrices can be calculated via 
the formula 


00 

E 

k=l 


ekz^ = 


(1-A?z2)(l-Aiz2) 


z-^S2 + ^"'51 + 


-li 


Multiplying by (1 — \\z^){l — X^z^) and equating the coefficients we obtain 


^2 =01 
=02 

So = 03 - (A? + A2)0i (21) 

0 = 04 — (A^ + A2)02 

0 = 0fc — (A 4 + A2)0A :—2 + A 4 A 20 A:—4) A: > 5. 


We note that also in the general case there exists a (more complicated) linear recursion for the 0^ 
matrices. Expanding cosh AA* in a Taylor series, combining with Proposition 1, we obtain 



4g^('^l+-^2) 

2^r 


k=l 


where Ck{u}) are trigonometric polynomials. Using the first few values of the coefficient functions 
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Cfc(w) and (21) 


Ci{uj) = (1 — cosw)^ 


C 2 {oj) = —isino; 


02 


C3(a;) = - 0i(A2 + Ai)-03 


2 + cos oj 


xsincj 

C4{uj) — -4^®4- 


Thus we may obtain a higher order expansion of the spectral density, e.g. 
/ma(‘^) = ^( ~ 2'S'2(1 - cosw) + A l"-2(Ai + A2)5'2(l - coscj) + isinw5i 


27r V 
— ( (1 — cos w) 


<50 + 'S'2(Af + A 2 ) 


+ S’2(Ai + A2)^ 


— 5o — i sin ui 


5i(Ai + A2)) +0(A3)^, 


from which the statement of Theorem 1 again follows. However, it is not clear how to obtain higher 
order expansions for the process itself. 


5 Proofs 


Proof of Proposition 1. Let r(t) denote the covariance matrix, i.e. T{t) = EYol^"''. Then 

T{t) = [ e‘*Vv(A)dA, t G M, 

Jr 


and applying Cauchy’s theorem componentwise we have for t > 0 

Tit) = ^J^e^^Riz)dz, 

where p C (—oo,0) x IM is a closed curve, which encircles the zeros of detP( 2 ;). Since T is a 
covariance matrix, r(— t) = r(t)'''. It is clear that the autocovariance function of the discrete 
process (TnA)neN is r(An), so by the inversion formula for discrete processes the spectral density 
can be calculated as 


^ 00 

= ^ E 


—ikio 


2t: 

1 

1 


r(AA:) 


k=—oo 
0 


g-ifca;r(_AA:)^ + Y e-'^‘^T{Ak) 


Lk=—oo 


k=l 


dvr^i 

1 

dvr^i 


n 00 n 00 

Jp ^^0 ’ 


1 


P k=l 
Az—iio 


( 22 ) 


LJp 


_ gA^+io; 




, wG(-7r,7r), 
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where the change of the sum and integration is justified, since Dlez < 0 on the curve p. 

Consider the Laurent-series 

1 _ 

1 _ = I] and 4(a;)z^ 

k=-l k=-l 

Note that either function has a pole at 0 only if a; = 0, therefore for cj 7 ^ 0 the series above are 
usual Taylor-series. In the following we assume that w 7 ^ 0. Adding the two expressions 

1 — e~^ sinh z 

1 — 1 — + e~^ — 2 cos oj cosh z — cos cj ’ 

which is an odd function of z, therefore in the series expansion 

C2fc4) +4fc(w) = 0. 


On the other hand 


+ e ^ — 2 e 


cosh z — e 


\ _ qZ iw Y _ g^+io; + e ^ — 2 cos LO cosh z — cos OJ 

is an even function of z, therefore 

d-2k+l{oj) - C2fc+l(w) = 0. 

Summarizing, we obtain 


= -1 -i 


sm w 


cosh z — cos OJ 


dk{oj) = 


For the coefficient Ck{oj) we have 


1 


sinh z 


2 cosh z — cos OJ 


y~lc2A;+lO 


OJ z 


2k+l 


(23) 

(24) 


A:=0 


and 


1 ^ i 


sincz 


2 cosh z — cos OJ 


^C 2 fc(a;)z^^, 


k=0 


so for k>l the coefficient C 2 fc(a;) is purely imaginary. From (24), using the notation of [9] formula 
(18) we see that C 2 fc+i(cz) = —Ck{oj)/2. 

Inserting the series expansion into (22) and using that R{z)'^ = R{—z) we obtain 


/a(w) = ^;^^A^ ^Cfc(w) J z^R{-z)dz + dk{oj) J z^i?(z)dz 

Then, changing the variables and using (23) and (5) we have 

Ck{oj) f z^R{—z)dz + dk{uj) f z^R{z)dz 
J p J p 

= {—l)^^^Ck{oj) f z^i?(z)dz + (—l)^^^Cfc(a;) f z^R{z)dz 
J —0 J 0 


(25) 


( 26 ) 


= (- 1 )^+' 4(4 / 

J pU—p 

= (-l)^+^Cfc(w)27ri0fc. 


z^R{z)dz 
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Substituting into (25) 


/a(w) = — ^(-A)*'0fcCfc(a;). 

k=0 

Taking into account that 0^ = 0 for A: < 2(p — g) — 2, the proof is ready. 
Proof of Proposition 2. We have 


R{z) = 


adiP{z)Q{z)T,L{&diP{-z)Q{-z)) 

detP{z)detP{—z) 


detP{z)detP{—z) 


Since the degree of the numerator is less than that of the denominator, for the partional fraction 
decomposition (11) holds for some matrices a{£),(3{£) G M^. Note that 5(—z)"*" = S{z) implies 
/3{£)'^ = a{£). By simple properties of the residue the second summand in (22) is 


r gAz-icj ''“If 

= A / I 

p /?_1 p 


27ri /n 1 — 


( a{£) m \ 

_eA.-i^ V^-A^ + -z-Xe) 


pAXi-iu) 

X/ 1 _ gAA^-ia;“(^)- 


In the same way 


Therefore, by (22) 


27ri Ll- 


pd 


1 P‘^ / AA^-iw 1 \ 

A(-) = ^ E ( 1 _ + i_eAA, . 

i=\ ^ ^ 

from which simple manipulation gives (10). □ 

Proof of Proposition 3. In the general case the partional fraction decomposition of R{z) reads as 
(13), with some matrices a{£, j), j3{£, j) € M^. Similarly, as in the previous case we obtain 


27ri .L 1 — 


-i?(z)dz = 


f=i j=i 


pAz luj \ a{£,j) 


where D stands for differentiation. Noting that 


pAz-XiV ^ ^ I 

_ pAz—iuj ^ _ pAz—iuj ' 


one can show that for j >2 


2)0-1). 


1 \ _ sf{u},Xe) 

1 — )z=\ (cosh AAf — cos cj)-i ’ 
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with s^{io,X£) being a trigonometric polynomial in u of degree j — 1, whose coefficients depend 
on A and A^. Similarly, for j >2 




1 

\ _ qAz+xui 



sf{-u},Xe) 
(cosh AA^ — cos ujy 


Substituting back into (22) we obtain (12). □ 

Proof of Lemma 1. Recall definition (9). To ease the notation put h{z) = 1/(1 — e^"’"''^), and 
y = . Then, for the first few derivatives (all the derivatives are meant with respect to z) 


h'{z) 


y 

(i-y)2’ 


h"{z) 


y'^ + y Ufa/X _y^+ ^y‘^+ y 

(1-2/)3’ (l-y)4 


In general, using induction it is easy to see that 


h^^\z) 


y^niv) 

(l-y)n+l’ 


n = 1,2,.. 


where An is a polynomial of degree n — 1, for which the recursion 


An+i{y) = {y- y‘^)A'n{y) + An{y){ny + 1) 


(27) 


holds. These are called Eulerian polynomials (not to be confused with Euler-polynomials). The 
coefficients are the Eulerian numbers, i.e. An{y) = A{n,n — l)y^~^ + A{n,n — 2)y”“^ + ... + 
A{n,0). The combinatorial interpretation of the Eulerian numbers is that A(n,k) is the number 
of permutations of {1, 2,... , n} in which exactly k elements are greater than the previous element. 
From (27) we obtain 


A{n + l,k) = {k + l)^(n, k) + {n + 1 — k)A{n, k — 1). 

Induction gives that A{n,n — 1) = A{n,0) = 1, and 

A{n, k) = R(n, n — 1 — k), /c = 0,1,... , n — 1, (28) 


that is An is a symmetric polynomial. 

Since (1 — e'‘^)(l — e~^) = 2(1 — coscu), from (9) we have 


n! Cn{oa) = h^"'\o) 


e^‘^An{Pn 
(1 - 


[2(1 — costi;)]"-+^ 


For odd indices, n = 2k — l,k = l, 2,..., using (28) 


A2k-i (e^) = A{2k - 1,+ A{2k - 1, + ... + A{2k - 1, l)e‘^ + yl(2A: -1,0) 

= [A{2k - 1, 0) cos(/c - l)a; + A{2k - 1,1) cos(/c - 2)w + ... + 2"^yl(2/c - 1, fe - 1)] . 
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The second factor is a polynomial of coscj of degree k — 1. For the first factor — e = 
—2(1 — cosw), therefore 


—\u)\2k 


{2k - l)!c2fc-i(a;) = 




[2(1 — costi;)]^*^ 

2(-l)*^ [A{2k - 1,0) cos{k - l)w + yl(2A: - 1,1) cos{k - 2)uj + ... + 2-^A{2k -l,k 

[2(1 — cos w)]^ 


(29) 


Forn = 2A:, A: = 1, 2,..(28) implies A2fc(-1) = 0, i.e. A 2 k{y) = (1 + y)^ 2 fc-i(l/), where 
A 2 k-i{y) = A{2k - 1,0)y2^-2 + I(2A: - 1, l)y2^-3 + ... + A{2k - 1, l)y + A{2k - 1, 0) 


is again a symmetric polynomial of degree 2k — 2. Thus, using the calculation above, and that 
(1 + e'‘^)(l — e~^^) = 2isina; we obtain 


{2k)'. C2k{u}) = 


e‘^A2fc(e‘"^)(l 

[2(1 — cos w)]2^+^ 


4(—l)^isin w yl(2A: — 1,0) cos(A: — l)u} + ^(2A: — 1,1) cos{k — 2)u! + ... + 2 ^A{2k — l,k — 1) 


[ 2(1 — 008 0 ;)]^+^ 


(30) 


Apart from constant factors the statement is proved. 

Expressing cos nui as a polynomial of cos oj serves as a definition of the Chebishev polynomials 
Tn, i.e. 

cosnu = Tn{cosuj) = ^ ^(-1)^ ^^^_^^^,' (2cosa;)^-^^ 

From this we see that the coefficient of (coscu)*^ equals 2”“^. Thus the coefficient of (coscu)^”^ on 
the right-hand side of (29) is (—1)^2^“^, from which we obtain the value of the leading coefficient. 
After noting that A{2k — 1,0) = A{2k, 0) = 1, the value of the leading coefficients follows similarly 
in the even case. Finally, from (27) we see that An{l) = n!, from which the formula for the product 
of the roots follows. 

Thus we have shown that the polynomials qk-i,rk-i defined via 

qk-i{cosuj) = {-l)^{2k - 1)![2(1 - cosa;)]^C 2 fc-i(a;), 
isina;rfc-i(cosa;) = (—1)^(2A:)![2(1 — cosa;)]^'''^C 2 fc(a;), 


have the stated properties. From (29) we see that qk-i and Vk-i are linear combinations of 
Chebishev polynomials 

qk-i{x) = 2 [A{2k - 1, 0)Tk-i{x) + A{2k - 1, l)Tk- 2 {x) + • • • + 2-^A{2k - 1, A: - 1)] , 
rk-i{x) = 4 \A{2k - 1,0)rfc_i(x) -f A{2k - 1, l)Tk-2{x) + ... + 2~^A{2k - l,k -1) 
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Now we turn to the statement about the roots. Indirectly assume that qt-i has a real root 
in (—1,1), say coscuo- Then from (29) we see that 742 ^- 1 ( 6 “^°) = 0. This is a contradiction, since 
Frobenius showed in 1910 that the roots of the Eulerian polynomials are real (for a recent work on 
roots of generalized Eulerian polynomials see Savage and Visontai [20]). Similar reasoning shows 
that rfc-i(x) has no real root in (—1,1). □ 


Proof of Theorem 1. 
As 


The proof relies on analyzing the spectral density /ma(‘^) process X^. 


00 

cosh AAj — cos oj = 1 — cos oj + 

t=i 


(2£)l ’ 


using Corollary 1 and (15) the asymptotics of the spectral density of the moving average process 
is 


i=i 


■ypd 


pd 


2tt 

/\2(p-g)-l 


— COSCJ 




k=0 


2tt 


2P‘^{1 - COSUj)P^C2(p-g)-l{uj)Q2ip-q)-l 


as A 4 , 0. Write 


(31) 


A2{p—q)—l 

/ma(‘^)- ^ -[2(1 - cosw)]^"*“^^“'') [2(1 - cosw)]P"‘?C2(p_g)_i(a;)02(p_q)_i. (32) 

It is clear that in (32) the factor [2(1 — coscj)]^(‘^“^)+'^ corresponds to {Id — in the 

moving average representation. 

For the factorization of the remaining terms we need that 


{1 - r)P‘^){l - r]e-^‘^) 


2r] 


^1 — COSCJ + 


2r, ) 


thus solving the equation = (1 — r/)^/(2r/) we have for the solution 


r/(0 := 1 - C ± - 2?, 


where the sign is chosen so that \'r]{f)\ < 1. This is possible, since the product of the two roots is 
1. That is 

1 - coscj - C = ^^^(1 - ^(Oe'‘^)(l - V{^)e~'‘^)- 
Therefore combining the above with Lemma 1 we obtain 


[2(l-cosa;)f ''c 2 (p_g)_i(w) 


^p-q-l IfjJl \l-COSU-C2{p-q)-l,j) 

{2{p -q)- 1)! 

nj=r^(l - Vi^2(p-q)-l,j)e'‘^){'^ - V{f.2ip-q)-l,j)e~'^) 
{2{p-q)-l)inW^v{C2ip-q)-i,q) 
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We conjecture that the zeros ^ 2 k-i,j are all real and greater than 2. This is true for k = 
1,2,..., 8, however we cannot prove it in general. For real zeros the r]^s are also real (we did prove 
that ^ 2 k-i,j ^ (0)2)), thus in the factorization everything is real. When there is a non-real root ^ 
then necessarily its conjugate ^ is also a root, and one can check easily that rj{^) = ry(^), therefore 
in the factorization the coefficients are real. □ 
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